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ABSTRACT 


The diffusivity matrix [Dj in the ternary liquid 
system comprising ^ polystyrene, toluene and cyclohexane 
has been determined employing the diaphragm cell method. 

The evaluation of the diffusivities from the mass concen- 
tration vs. time data was done by a numerical multi- 
variable optimization technique, namely, the modified Box 
method. Stronger multicomponent effects as well as 
concentration dependence of the diffusivities were observed 
compared to the data r ported in the literature for the 
system with polystyrene of higher molecular weight. 
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Chapter 1 
INTRODUCTION 


Molecular diffusion arises due to the movement of 
individual molecules by virtue of their thermal energy. 
Diffusion phenomenon in fi^ases is described by the kinetic 
theory. However, the diffusion in liquids and solids are 
more complex in nature than gases and are not well under- 
stood. 

The phenomenon of molecular diffusion ultimately 
leads to a completely uniform concentration of species 
throughout a mixture, which may initially have been non- 
uniform, The diffusive mass flux of component A in a 
binary mixture is given by Fick’ s first law as: 

h = -®AB (I-') 

where is the mass concentration gradient for the species 

A. Eqn. 1.1 gives the definition of binary diffusivity 

Multicomponent diffusion is characterised by the fact 
that the flux of one componoit is influenced by the 
concentration gradients of other species except* for ideal 
cases. Pick’s first law can be generalised for a multi- 
component mixture as 

( 3 ) = - 


[Dj ( vP ) 


( 1 . 2 ) 
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where, for a mixture containinc^ n ccmponents, ( 3 ) and 
(vP ) are n.~1 dimensional column vectors of the mass flux 
and mass concentration scradient^ respectively and [dJ is 
an n -1 dimensional square matrix of diffusivities. 

The components of the matrix [dJ , i.e,, the 
are termed as multicomponent diffusion coefficients and 
are different from the binary diffusivity in Eqn.l , 1 . 

The relation between the fluxes and the mass concen- 
tration gradients in a multicomponent mixture is also given 
by the generalized Stefan-Maxwell equation [Cussler and 
light foot, l965aJ: 



In a multicomponent system, it is possible to observe 
various “interaction phenom^a” [Toor, 1957; 1964], such as 

(i) osmotic diffusion (transfer of a component in the 
absence of composition gradient for that species), 

(ii) diffusion barrier (no transfer of a component even 
though a compos itipn gradient exists for it) and (iii) 
reverse diffusion (diffusion of a sp ecies against its 
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composition e:radient). These different interaction phenomena 
are attributed to what are grossly termed as "multicomponent 
effects". 

Most of the chemical engineering processes, for 
instance, distillation, absorption, condensation, evaporation 
etc, involve multicomponent mass transfer; and diffusivity 
data are required for their analysis and design. Most 
often one finds that the experimental data are lacking, 

Diffusivities of polymer solutions are required for 
the analysis of many polymer formation and fabrication 
processes. The rate of polymerisation is strongly influenced i 
by the molecular diffusion of monomers, initiators, long 
chain free radicals, dead polymer chains or low molecular j 

weight condensation products. The processes such as 
separation, purification, introduction of plasticisers and 
other additives, making of foams etc. also involve mass 
transfer of polymer molecules in different medium and the 
diffusivity data are required for the analysis of these 
processes as well. Unfortunately, data in the literature 
are virtually non-existent for ternary and multicomponent 
systems involving polymers. 

The aim of the present work was to obtain Wiperimental 
diffusivity data in ternary polymer-solvent system and 
to examine the multicomponent effects. The Diaphragm cell ' 
method [Cussler, 1976 ] was employed for the measurements. 



The system was polystyrene-toluene-cycloh(S£ane, which 
was chosen for its strong multicomponent behaviour, as 
reported by Cussler and ligbtfoot [l965aj , 
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Chapter 2 
LITERATURE EEVIBW 

In this chapter, a brief discussion is made on the 
various techniques for experimental determination of 
binary and ternary diffusivities. This is followed by 
discussions on the various experimental data reported in the 
literature and the theories for prediction of diffusivities 
in polymer-solvent systems. 

There are a wide rahe:e of methods for measurement of 
diffusivities in liquids. The use of different types of 
equipments and the« methods have been discussed by 
Oussler ] . More extensive discussion can be found in 

the book by Dunlop et al [l972j . The mathematical analysis 
to obtain diffusivities from the raw data are discussed by 
Crank [l 975 j. 

The more frequently used methods employ either the 
-a 

Cojdy interferometers or the diaphraepi cells of different 
constructions. The Gouy interferometer [Costins; et al, 1949] 
measures the refractive index gradient between two 
interdiffusing solutions with respect to position and time. 
The diffusion coefficients derived from Copy interferometric 
studies are highly accurate, often better than 0.1% 

[Cussler, 1976J . Most of the studies on polymer- solvent 
systems, reported in literature, have been done with the 
Gouy interferometers. 
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In the diaphr^m cell method, the variation of concen- 
tration of two intend iff using solutions separated by a 
porous diaphragm in the cell is measured with time. These 
data are utilized to determine the diffusivities as will 
be discussed in the Chapters* 


The effective diameter of the pores in the diaphra.crm 
must be such that gross streaming through the diaphragm 
is avoided, but, at the same time, the pores must be large 
in comparison with molecular dimensions so that the diffusion 
will take place under conditions comparable with those 
obtained in free diffusion. Diaphragms made of 2-3 mm thick 
fritted glass had been used to determine diffusivity of 
KOI in aqueous solution [Gordon, 1945 ] . However, thin 
diaphragms made of different synthetic materials are used 
now-a-days, which are more suitable for using with organic 
solvents and polymers. The design of diffusion cells, 
since its earliest uses in the forties, has also been 
modified according to requirements. A later version of a 
diaphragm cell has been described by Robinson and Stokes 
[I 959 J and a more recent version was used by Geankoplis et al 
[1978] to measure diffusivity of urea and KOI in albumin 
solution, 

A comparison of the Gouy interferometer and the 
diaphragm cell [Cussler and Dunlop, 1966] asserts that no 
theoretical advantage in either of them should be expected. 
The mathematics involved in both the cases has got similar 
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de^ee of difficulty. The principal advantage of the 
diaphrafifoi cell is its simplicity of construction and use. 
Therefore, even though the Gro^y interferometer is often 
more accurate, the diaphragm cells are used when ease of 
obtaining data is one of the concern .5, 

Thersf; are two techniques for the analysis of raw 
ternary data to determine diffusion coefficients [Cussler, 
1976 J. One is the pseudo-binary, method [Fujita and Costing, 
1960J Eevzin, 1972 J, used commonly with the Ooidy inter- 
ferometers. 

The other technique based on numerical minimization 
of a suitably foimulated objective function, is discussed 
in detailn in Chapter 4. A numerical technique, based on 
the non-linear least-squares is described by Patel et al, 

[1969 J. 

Experimental data for polymer- solvent systems are 
extremely meagre. One the earliest experimental work 
was reported by Tsvetkov and Klenin [I958j» who studied 
the diffusion of two grades of polystyrene (M.^ = 3,500,000 
and 5,000,000) in different solvents such as butanone and 
car bontetra chloride. The differential diffusion coefficients 
were determined by keeping the concentration diff erence as 
small as 10 ^g/ml in polarisation interf eromater. 

.The concentration dependence of diffusivity, 
described as 

D = Dq (1 + kjjP) (2.1) 
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where p is the concentration of polystyrene in g/ml and 
and kjj are constants, was also studied. The findins:s are 
reported in Table 2.1, 

Gussler and Lia:htfoot [1,965a] reported binary as well 
as ternary diffusivity data for polystyrene (M^ = 193,000) 
in toluene and cyclohexane. Measurements were done by a 
G-0% interferometer at a temperature of 28.00°C, Actual 
diffusivity values for different compositions were found 
by repeating all experiments at the same average concen- 
tration but different concentration differences and extra- 
polating the results to zero concentration difference. Both 
the binary and the ternary diffusivitias showed strong concen 
tration dependence. The ternary coefficients showed large 
deviations from Pick's law. The cross term coefficients 
(1)^2 ^21^ were found to be much lars:er than the main 

term coefficient and ^22^’ seen ftom the 

data reported in Table 2.2 and Table 2,5. 

Ternary diffusion of mixtures of two monodisperse 
polystyrene of different molecular weights = 193,000 
and = 234,000) in toluene was studied by Gussler and 
Lightfoot [1965b]. The results also showed large multi- 
component effects. The results were in qualitative agreement 
with those values predicted from the extensions of hydro- 
dynamic models [Groot and Mazur, 1962; Scott, 1949 | Plory 
and Krigbaum, 1950 j , The results of the experiments 
are reported in Table 2.4. 
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Table 2,1 ; Concentration dependence of binary 

diffusivity in polystyrene- solvent systOE 
[Tsvetkov and Klanin,1958 j 

Note; and refer to Sq.n,2,1 


Mol.wt, 

I 

{ Solvent 
! 
t 

t 

— f 

{ L- 

t 2 ® 

I (cm /sec.) 

r 

^T) 

I cm^/g.) 

3.5x10^ 

But anon e 

1 .Ixio"^ 

1 .6x10^ 

3.5x10^ 

COl^ 

0.45x10"^ 

1 . 5x1 0^ 

5.0x1 0^ 

CCl^ 

0.33x10"'^ 

6.0x10^ 


Table 2,2 : Binary diffusion coefficient of polystyrene- 

solv^t system 

[Cussler and Light foot, 1965a] 


Mol.wt 

f 

t 

{ Mass fraction 

n 

f 

J Solvent 

; Bif f usivityjB^g 

1 (cm /sec) 

I polymer 

f - - 

1 

f 

t- - - 

193,000 

0.05 

Toluene 

8.752x10“^ 



Cyclohexane 

I.OIOxlo”^ 
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Table 2.3 : 

Ternary diffusivity of polystyrene-toluene' 
cyclohexane mixture 

LCusaler and Light foot, 1965a 1 

Set 1 


Components ; 

1 = Polystyrene (M = 193, §00) 

w 

2 = Cyclohexane 


3 = Toluene 


T emary mixture 

1 

2 

q 


0.0500 

0.0500 

¥ 


0.0500 

0.9000 

2 


9000 


3 


O.O 9 QO 

e .0500 

^x1 0 ^ 

(cm^ / sec) 

8.9 

10.0 

3 ^ 2 X^ 0 ^ 

{ca? / sec) 

-1 .6 

1.3 

Lgi^io'^ 

(cm^/sec) 

- 8.9 

169.6 

D 22 XIO'' 

(cm^ /sec) 

205.1 

202 .0 


Continued . 
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Table 2.3 (Continued) ; 

Set 2 

Components ; 1 = Polystyrene (M^ = 193,000) 

2 = Toluene 

3 = Cyclohexane 


Ternary mixture 

1 

2 


0.0500 

0.0500 

Wp 

0.4751 

. 0.0503 


0.4749 

0.8997 

7 2 

^xIO (cm /sec) 

5.8 

1.8 

7 2 

Di 2 '’x 10 (cm /sec) 

-6.7 

-6.7 

I>2i3c 10' (cm /sec) 

14.5 

-0.9 

(cm^/sec) 

171.8 

178.9 
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Table 2,4 ; Ternary diffusivity of polystyrene 1- 
polystyrene 2 -toluene mixture 

[Cussler and Lightfoot, 1965b ] 
Components : 1 = Polystyrene (M_ = 193,000) 

w 

2 = Polystyrene (M^ = 234,000) 

3 = Toluene 



■n 1 

|vi+W2 1 
! 1 

V f 

t 1 

! 

D^^xio'^ i 
(cm^/sec) } 

f 

r 

i 

(cm^/sec) 1 

y. 

Dg^xio”^ 
(cm^ / sec){ 

(cm /sec) 

E>:p eri- 

0,0000 

3.77 

0.00 

0.00 

3.45 

mental 

0,0230 

4.46 

1 .14 

1 .73 

4.69 


0.0500 

4.20 

4.17 

4.21 

4.23 

Theore- 

0.0000 

3.77 

0.00 

0.00 

3^45 

tical 

0.0250 

4.43 

2.94 

2.72 

4.32 


0,0500 

5.10 

5.92 

5.47 

5.32 1 
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Binary diffusivity data for polystyrene in cyclohexane 
from the Gooy interferometric studies were reported by 
Cowie and Cussler [196?], which showed variation of diffu~ 
sivity with molecular weight of polystyrene. The measured 
diffusivities were actually those averaged over the poly- 
dispersity of the polymer, which was, however, very low. 

Results are reported in Table 2.5. 

Diffusion of polyst^ene in cyclohexane at the dilute 
and the infinitely dilute range were studied by Yrentas 
and Duda [1976a,b J . Diffusivity data taken from various 
sources were examined against those predicted on the basis 
of the Kirkwood -Riseman theory [Kirkwood and Riseman, 194-8 
} lamakawa, 1971 J . The predictions and the experimental 
data were found to be in good agreement for polystyrene -■ 
cyclohexane system. 

Diffusion of ethylbenzene in molten polystyrene 
(M^ = 275jCOO) was studied systematically by Duda and 
Ni [1978J at the temperature range from 115.5 to 178°G 
and for concentrations as high as 70 weight percent ethyl- 
benzene. Experiments were carried out using a quartz-spring 
sorption apparatus [Duda et al, 1973j, which was modified 
for high tanperature measurements. The data show strong 
concentration dependence of diffusivity, which is more 
pronounced at low tanperatures and lower solvent concen- 
trations, Erom the comparison of these data with those 
predicted from the free volume theory [Duda and Yrentas 1977a,bJ 
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Table 2,5 : Binary diffusivity of polystyrene -cyclohexane 
mixture - variation with molecular weight of 
polymer 

Coussler and Cowie, 1967 ] 


Temp. 

(°c) 

"T" — 

\ Molecular 

{ weight 

! 

1 

-..1 

1""' 

j Diffusivity, Dj^^xio"^ 

1 (cm /sec) 

f 

1 

35 

1 94000 

3.082 


414000 

2.225 


929000, 

1 .436 

25 

181000 

3.494 


481000 

1 .816 


1 040000 

1.157 
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it is found that the theory consistently overestimated 
the diffusivities. 

The mass transfer in the concentrated poljnner 
solutions can not be described by the classical theories 
of diffusion. This non-classical behaviour is supposed to 
be associated with the sluggish relaxation of the large 
polymer molecules from the bulk. These behaviours were 
studied by Alfrey [19^5] . Work was later extended by 
Hopfenberg and Frisch [I969J and Duda [1975j. A diffusion 
Deborah number was introduced as 

( 2 . 2 ) 

where K is the characteristic time of the fluid and 
m 

is the characteristic diffusion time. The value of 
Deborah number for a given set of conditions identifies 
the type of diffusional transport that can be expected. 

When Deborah number is small (say less than 0,1), the local 
relaxation in the system is fast compared to the diffusion 
process. On the other hand, a large Deborah number (say 
greater than 10) signifies that there is essentially no 
time variation of the structure of the polymer during 
diffusion process. This type of diffusional transport was 
denoted as elastic diffusion [Vrentas and Duda, 1979 J . 

The characteristics of polymer solutions change 
significantly with the change in concentration . In the 
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infinitely dilute polymer solution, the polymer molecules 
are widely dispersed and there are no interactions between 
the individual polymer chains • As the polymer, concen- 
tration is increased, the polymer molecules begin to 
interact hydrodynamically with each other, even though the 
domains of the individual polymer molecules do not yet overlap 
This concentration region is referred to as the dilute 
solution region. At a further increased polymer concen- 
tration, the domains of the polymer molecules begin to overlap 
so that there is a considerable amount of polymer-polymer 
contact, and entanglements are formed between polymer 
molecules. This concentration region is denoted as 
concentrated polymer solution [Vrentas and Duda, 1979j« 

Although a general theory which can describe molecular 
diffusion throughout the complete concentration range is 
yet to be developed, independent theories have been 
proposed which are applicable to specific region of the 
concentration intervals. A brief discussion is given 
below. 

For dilute polymer solutions, the concentration 
dependence of mutual diffusion coefficients can be 
expressed as [Vrentas and IXida, 1979 j: 

= Do (1 ^ 2 ^ (2 *3 ) 


D 
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where p g is the mass concentration of the polymer in the 
solution, is the mutual diffusion coefficient at zero 
polymer concentration and kjj is a coefficient. Thus, 
in dilute polymer solutions, prediction of diffusivity 
reduces to determination of and at a particular set 
of conditions for the particular system. The determination 
of involves consideration of intramolecular forces, whereas 
both intra- and intermolecular forces are to be considered 
in the determination of Various methods for these along 

with references are discussed by Duda and Vrentas [1979 J • 

Molecular diffusion in concentrated polymer solution 
is a more complex phenomenon and is not well understood. 
Nevertheless, the free volume theory described by Fujita 
[1968] , which was further modified by Vrentas and Duda 
[1977a;, b ], rendered positive contribution in understanding 
the phenomena and also in prediction of diffusion coefficients. 

The predictive methods described in the precjfeeding 
section are based on approximate theories of molecular 
diffusion. To evaluate these methods and to implement 
necessary improvements, reliable diffusivity data are 
required. How ever, very little experimental work has been 
done to determine diffusivities in polymer- solvent systems, 
most of which is due, to Cussler and lightfoot [ r965a,b ]. 
Various notable experimental works on determination of 
diffusivities are listed by Cussler [1976 j . 



Chapter 3 

EXPERIMENTAL DES^AILS 


3.1 Theory of Diaphragm Cell 


The diaphragm cell used for the determination of 
diffuaivity is shown schematically in Figure 3.1. Lt 
consists of two compartments A and B, with provisions for 
stirring, separated hy a diaphragm, which can either be a porous 
plug or a membrane [Stokes, 1950^ 1951 j • Ttie compart- 
ments are filled with solutions of different concentrations 
and are kept well-stirred. These concentrations change 
gradually with time as a result of diffusion across the 
diaphragm. The change in concentration is measured at 

^ . a- 

different intervals of time and the data^-is used to compute 
the diffusion coefficient (s) as described below. 


For a binary solution, the mass balance equation for 
the solute within the porous plug or membrane is [Oussler, 


1976 J 






( 3 . 1 ) 


and the boundary conditions are 


x>0 , = ^lA 

X = 1 , P 1 * P-IB 


(3.2) 

(3.3) 


where 1 is the thickness of the diaphragm and ^ and 
are the mass concentrations of component 1 in compartments 



Well - stirred 
solutions 



FIG.3J 


DIAPHRAGM CELL DIFFUSION 
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A and B. 


When the volume of the pores in the diaphragm is much 
less than the volume of the entire diaphragm cell, a steady- 
state concentration gradient is very quickly established 
in the diaphragm. However, this is actually a "pseudo-steady 
state" value, which does change very slowly with time 
[Cussler, 1976 J • With this pseudo-steady state approximation. 


lqn.3.1 becomes 
0 . = - 


dx ^AB 




(3.4) 


Intonation of this equation and substitution in the overall 
mass balance equation yields 


' pO = 

IB " 


(3.5) 


where the subscript 'o' indicates the initial concentration 
and P is the o^l constant. It is related to the cell 
parameters as 


P = T (^ + ^) (3.6) 

^ ^B \ 

where A is the total area of the pores in the diaphragm and 
and are the volumes of the two compartments of the 
ceibl* , 

Bor binary mixtures with known diffusion coefficient, 

Hqn. 3. 5 may be employed to d etermine the cell constant . 
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Once the cell constant is known, the same equation can be 
employed to determine diffusivity of unknown mixtures. 

The pseudo-steady state assumption has been thoroughly 
tested by Mills et al [1968J and it was reported that the 
above approximation leads to an error less than 3 percent, 
even when the pore volume is 10 percent of the total 
volume. Since the pore volume is commonly less than even 
1 percent of the cell volume, the approximation is indeed 
very good. Also, the time required for attainment of 
the p seudo— st eady state within the diaphragm, t was given 
V fGordon, 19453 

iL 

^AB 


to = 1.2 


(3.7) 


It should be ensured that the experimental run-time t is 


such that t > > t^. 


The extension of the binary diffusion equations for 
the diaphragm cell for the ternary diffusion yields the 
following equations [Toor, 19<''4a,b; Cussler, 1976j; 




( Dii -'^ 2)^^1 + ^ i 2^^2 


a <- O- 
‘'l ""2 


exp (-p0^ t) 


O O 
2 1 


(3.8) 



22 


APa = 

4 “ 

where , 

= I* +D 22 ^ Di2^21^’ 

‘’ 2=7 »11 +®22 “l2“21> 

the symbol A represents the difference between the upper 
and the lower compartments and the \-j'3 ternary 

diffusion coefficients. They can be evaluated from vs. t 
data. Since there are four unknowns (the and only 

two equations, it is necessary to perform at least two experi- 
ments to determine the four D^^'s. However, according to 
Cussler [1976 J, twice the minimum number of experiments is 
desirable to obtain accurate values of diffusivities, 

3.2 Difinhragm Oell 

The diaphragm cell used in the present work was 
similar to the one described by Geankoplis et al [1978 j. 

A scale drawing of the cell is given in Figure 3.2. 

The cylindrical shell and the caps of the cell were made 
of nylon, as nylon is unaffected by the organic solvents used. 


^ 21 ^ 1 ^ ^^22 ~ 2 
1 ■" 2 


exp (“pUj t) 


1 exp{-M2t) 


(3.9) 



3x6 MM THREAD AND NUT 



19x4 MM HOLES 
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Table i Parts~list for diaphragm cell 

[Pef. : Fig. 3. 2 j 


Part 


Description 


1 4 mm nylon screw used as cap 

2 Magnetic stirrer 

3 Brass fastener 

4 Steel capillary with 0,5 mm bore 

5 Teflon gasket 

6 Nylon shell 

7 Teflon diaphragm with 0,5 pm pores 

8 Steel plates with 19x4 mm holes 

9 Nylon column to support magnet 

^0 NiilffiAt 
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The diaphragm used was made of teflon and was obtained from 
m/s, Millipore Corporation, USA. According to the manufa- 
cturer’s specification, it had straight pores of 0.5 H m 
size. The diaphragm was housed between two stainless steel 
plates with coinciding large holes so that a large portion 
of the diaphragm was exposed to the solutions both the 
sides. Two cylindrical teflon-coated magnetic stirrers were 
placed in the two compartments of the cell, A nylon column 
was provided in the lower compartment to support the magnet 
near the diaphragm and a similar column was provided in the 
upper compartment to' maintain geometric similarity. The 
magnetic stirrers were operated by rotating two powerful 
horse -shoe magnets suspended from a frame outside the cell. 
Screw -threaded nylon caps were used to plug the holes in 
the compartments, those were used for transfer of solutions 
in and out of the cell. A narrow stainless steel capillary 
tube was provided in the upper compartment to allow overflow 
of the excess liquid while the hole is being plugged. A 
drawing of the set-up is given in Pigure 3.3. 

3, 3' Ixnerimental Procedure 

The two compartments of the cell were filled with 
two solutions of different composition as follows. The 
denser among the two was introduced in the bottom compart- 
ment to prevent mixing due to gravity. 



1. stirrer motor. 2. Temp, controlled water bath 
3. Horse- shoe magnets- 4- Diaphragm cell. 

5. Base, to which the celt is fastened. 


FIG. 3.3 SCHEMATIC SET-UP 
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First, the cell was inverted and the denser solution 
was introduced into the bottom compartment (which was 
actually at the top then) with a syringe. After filling, as 
the nylon stopper was tightened, the excess solution was 
passed through the porous diaphragm to the other compartment. 
The absence of air bubble was checked by visual observation 
against strong light. The cell wall was actually translucent 
and air bubbles were visible if they were present. 

Now the other compartment was rinsed quickly and then 
filled with the lighter solution. After filling, the nylon 
stopper was inserted and as it was tightened, the excess 
solution overflowed through the capillary tube provided 
for this purpose. Immediately after the capillary tube was 
covered with a teflon cap. 

Than the cell was placed in a constant temperature 
’Mater bath maintained at 28 + 0.5°G and stirring was started. 
It was found that the magnetic stirrers inside failed to 
rotate at high rpm. The stirrer was operated at 80 rpm 
for all the runs. After a predetermined time, the cell 
was taken out of the bath, the solutions from the two 
compartments were transferred to sample bottles using a 
syringe, and were subjected to subsequent analysis to 
determine their composition. This is described in the 
following sections. 
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5*4 Determination of Cell Constant 

E^n..3.5 was employed to determine the cell constant 
using binary mixtures of toluene and cyclohexane, the 
diffusivity of which was reported by Gussler and Lightfoot 
[ 1965a]. 

A calibration curve of mass fraction of toluene (w^) 
vs. density (P ) of toluene -cyclohexane binary mixture at 25°G 
was found. The densities were measured with a 20 ml 
specific gravity bottle using a Mettler H20 balance, which 
can weight with accuracy upto 5 decimal places of a gram. 

The calibration and the curve are given in Appendix A. 

The calibration curve was fitted to the polynomial 

Wj = P ^ (5.12) 

The coefficients were determined by non-linear least -square 
technique and are given belowt 

= -1 .8681 
* 6.2674 

« -7.0672 

a^ = 2.6948 (3.13) 

The initial and the final mass fractions of components 
in the solution were detemined using Bqn.3.12 from 
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dexisity measurements. Mass concentrations were consequently 
determined. The cell constant was calculated from three 
runs with toluene-cyclohexane mixture and the average value 
was used for the evaluation of ternary diffusivities. The 
determination of cell constant was done every time the 
diaphragm was changed, 

3*5 Determination of Composition of Ternary Mixtures 

It is possible to determine the diffusivities of a 
ternary mixture from the measurement of initial and final 
densities of the mixtures* But this method requires 4 sets 
of experimental data to obtain the diffusivities [ Cussler 
and Dunlop, 1966 j . However, if the mass concentrations of 
the initial and final solutions are known, then only 2 sets 
of data are required. Hence, in the present work, mass 
concentrations of the individual ^ecies were determined 
with the help of calibration as described below. 

Densities were measured for known composition of 
polystyrene-toluene-cyclohexane mixtures and a three- 
parameter calibration curve was made. The data and the curve 
are given in ilppendix A, The mass fraction of toluene w^ 
in the ternary mixture was fitted to the polynomial 

3 ' 3 

Wj = S S 3 .^ WpS ^ ^ 

i=0 3=0 


(3.14) 
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where Wpg is the polystyrene mass fraction and p is the 
mixture density at 25°G, The MG subroutine S{^2DAP, which 
employs the non-linear least-square technique, was employed 
to obtain the coefficients and are given below: 


1.11T'5 

a2.,= -3.8026 

^31- 

-18. 9424 

V'^-23.9294 

-12=58.3285 

a22=41.1453 

a^2= 

45.5083 

842=34.2428 

-,^=78.2732 

a23=79.6377 

a^^= 

67.7887 

a,, =62. 0261 
43 

..,4=114.1580 

824=105.4090 

^34= 

98.7856 

a^,j=93.4207 


(3.15) 

To analyse unknown samples, the density of the sample 
was measured first. Then a known amount of mixture (15 ml) 
was allowed to dry on a petri dish in an oven maintained at 
90*0. When all the volatile solvents were removed (this 
was confirmed from the constancy in weight at repeated 
weighing), the weight of the polymer residue was determined. 
Prom this, the mass fraction of the polymer was directly 
obtained. Prom the knowledge of P and "Wps? 
determined using Bqn.3.14 and the mass fraction of cyclo- 
hexane, W(,h, also determined. Mase ooneentratlons of 

all the components were subsequently calculated. 

The data as well as the calibration curve is given 
in Appendix A- The determination of diffusivity matrix [Pj 
from those data is discussed in the next chapter . 
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Chapter 4 

BESULTS AMD DISCUSSION 

In this chapter the d et emination of cell constant, 
verification of the reliability of the diaphragm cell and 
the ternary diffusivities are presented. 

4.1 Determination of Cell Constant 

Two diaphragms were employed in the present study, 
and the cell constants with these diaphragms were determined 
using a test solution of toluene and cyclohexane. The 
mixturedensity (P ) and mass concentration of toluene (P j) 
vs, run-time data are reported in Appendix B and the results 
are summarized in Table 4.1, 

From Table 4.1, it can be seen that the cell constants 
obtained were nearly the same with different durations of 
run (run-time) and also for the two diaphragms. 

4.2 Determination of Binary Diffusivities 

Binary diffusivity of benzene- cyclohexane and 
polystyrene -toluene system were determined and the values 
were compared with the values reported in the literature to 
check the reliability of the present method. The density/ 
mass concentration vs. run-time data are reported in 
Appendix B and the results are summarized in Tables 4.2 and 

4.3. ■ ■ ' : 
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Table 4«1 


Cell constants determined using toluene- 
cyclohexane binary mixture at 28 °C 


Average Composition 


Diffusivity 

[ Oussler and 
Light foot , 1 96 5a J 


Wij, = 0.5403 

W = 0.4597 
CH 

-5 2 

1.944x10 cm /sec. 


T 1 — -T 



{ Run-time { 

i (hr) j 

f 1 

- 1 - 1 

Cell constant 

(cm ) . 

Average cell 
constant 

(cm ^ ) 

Diaphragm 1 

0.5 

12.19 



1 .0 

12.50 

12.35 


1.5 

12.36 


Diaphragm 2 

0.5 

12.70 



1 .0 

12.12 

12.39 


1.5 

12.35 
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Table- 4. 2 ; Binary diffusivity of benzene -cyclohexane 
mixture at 28 °C 

Average composition ; Wg = 0.7875 

= 0.2125 

i.e., 20 mol % cyclohexane 
Cell constant = 12.55 cm”^ 


Run-time 

(hr) 


Diffusivity, 
(cm /sec. ) 


Average D 


AB 


( cm / sec . ) 


Reported D^ at 25 °C 
[Mills, 1 96 5j 

(cm^/sec, ) 


0.5 

1.944x10“^ 



1 .0 

2.185x10"^ 

2.03x10"^ 

1 . 90 x 10 "^ 

1.5 

1.956x10”^ 







Table 4.3 : Binary diffusivity of polystyrene-toluene 
mixture at 28 °C 

Average composition : = 0.05 

Wij = 0.95 
-2 

Cell constant = 12.35 cm 


Molecular weight, 


Diffusivity, D^ 
( cm^ / sec . ) 


32,000 


1 .039x10 


-6 
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The diffusiYity yalue for benzsne- cyclohexane system 
can be seen to be in good agreement with the value reported 
by Mills [1965] and the small difference may be attributed 
to the difference in temperatures (5°G) at which they 
were determined. 

The diffusivity for poly styrene- toluene mixture 

was found for polystyrene having a weight average molecular 

weight (B/^) 32,000, It is known that the diffusivity 

depends on the molecular weight of polymer. Data for the 

polymer of this particular molecular weight was not 

available for comparison. However, Vrentas and Duda [1979] 

gave the dependence of diffusivity for polystyrene and 

ethylbenzene mixture, for the polymer having molecular 

weight ranging from 1x10^ to 10x10^. Extrapolation of 

this data to the polymer molecular weight 32,000 gave a 
-6 2 / 

value 1x10 cm /secj which is in close agreement with the 
one obtained in this work. Since nature of toluene. as 
solvent is close to etl^ylbenzene, the diffusivities are 
not expected to be too much different. Moreover, Cussler 
and Iiightfoot [l965aj reported a value of diffusivity 
8.752x10’"'^ for polystyrene-toluene mixture of same 
composition for a polymer of molecular weight 193000. The 
diffusivity is approximately 85 percent of the value obtained 
in this work. This is again agreeable considering the 
large difference in the molecular weight, ^ ^ ^ ^ 
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Thus the cell and the method employed were considered 
satisfactory for the determination of diffusivities of 
binary and ternary mixtures, 

4,3 Determination of Ternary Diffusivities -Analysis of Pata 

To determine ternary diffusivities, experiments were 
performed with polystyrene-toluene-cyclohexane mixtures at 
three different compositions. Dor each of these compositions, 
4 to 5 independent experiments were performed with different 
run-times. The run'-times were chosen such that there was 
a significant canposition difference between the initial 
and the final values. Dor the sake of convenience, initial 
concentration differences ( P ^ = P ” P ±- q ^ between the 
upper and lower compartments were not varied from experiment 
to experiment, although this was not required. The samples 
were collected after each run and were analysed for their 
mass concentrations as described earlier [Section 3.5 L 

To determine the components of the diffusivity matrix 
from P^ vs t data, the following objective function was 
chosen; 

n 2 2 

1=2 S [ AP. (expt) -hP (calc)J (4.1) 

i=1 ^ j 

where hP^ (expt) are the differences in the mass concen- 
trations of component i between the upper and the lower 




37 


conipai“'tnieiibs of the cell at time t, A P ^ (calc) are those 
calculated using E(3.n.3.8 and 3.9 with a guess matrix [dJ and 
n is the number of experiments. 

The objective function was subjected to minimisation 
to obtain the matrix [dJ* The minimization was carried out 
using a subroutine (CPX) developed by Hughes and Malik [l979i . 
It utilised a multivariable search technique introduced by 
Box [Beveridge and Schechter, 1970 J . Domains of the D. .'s 

i tl 

were supplied in the foim of upper and lower limits. The 
program chose 7 different combinations of D. . values at 
random within these domains. The worst combination (for which 
the objective function is the largest) was replaced by a new 
combination selected by a specific criteria. The process 
continued till the minima for the objective function was 
determined with the desired accuracy, which can be set by 
multiplying the objective function by a suitable large 
number. The program which utilises the subroutine (GPX) 
is listed in Appendix C. 

The proper choice of the upper and the lower limits 
of the required a trial and error procedure, Pirst, 

the domains were kept wide and some rough minima was found. 

The minima was refined by narrowing the domains gradually 
such that in course of minimisation, the test D^^ values 
never hit the boundaries. 
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However, the nature of the function shows that there 
can he multiple minima. The local minima correspondirg 
to the feasible region was actually found. Thus a widd.y 
different set of D. . values might also minimize the objective 

X J 

function but the D. . values which were realistic has been 

X J 

found out using this procedure, 

4.4 Ternary Diffusivities-Results and Discussion 

Ternary diffusivity values were obtained for poly- 

styrme (M = 32,00o) - toluene- cyclohexane mijcture at 
w 

constant polymer concentration (Wpg = 0.i5) but three 
different toluene and cyclohexane compositions [Table 4.4]. 

Determination of ternary diffusivities require one 
of the components to be designated as "solvent” or the 
"third component". In the present work, there were two 
choices: first, cyclohexane and second, toluene. Two sets 
of diffusivities (Results - Sets A and B) were determined 
considering each of them as the solvent. 

It can be seen from Table 4.4 that the deviations from 
lick's law are large. The multicomponent effect, i.e., 
the dominance of the cross terms (D^p ^21^ 
main terms (D^ ^ and general stronger than those 

found by Gussler and Dightfoot [l965a j [Chapter 2] for 

polystyrene of molecular weight 193 jOOO» 



Tabl© 4.4 : Multicomponent diffusiTities of polystyrene— 
toluene-cyclohexane mixture at 28 °G 

Re suit -Set A 

Components ; 1 = Polystyrene (M = 32000) 

w 

2 = Toluene 

3 = Cyclohexane 


Ternary mricture 

1 

2 

3 



0.05 

0.05 

0.05 

^2 


0.90 

0.70 

0.475 

w 


0.05 

0.25 

0.475 


(cm^/sec) 

3.032x1 0 ”^ 

4.888x1 o“® 

3. 568x1 0 "^ 

1 1 

2 

(cm^/sec) 

1 .184x10-^ 

-2 . 1 53x1 0 ”^ 

-1.371x10”^ 

“21 

(cm^/sec) 

4.843x10 

-6 

9 . 992 x 10 ”^ 

-6 

4. 543x1 0 “^ 

—6 

“22 

(cm /sec) 

1 .719x10 

4.685x10 

4.499x10 


Continued. ... 
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T abl e 4 . 4 ( continued : 

Result- Set B 


Gomponeuts : 1 = 

2 = 

3 = 

Polystyrene 

Cyclohexane 

Toluene 

(M^ = 32000 ) 


Ternary mixture 

1 

2 

3 


0.05 

0.05 

0.05 

W2 

0.05 

0.25 

0.475 


0.90 

0.70 

0.475 

2 

^ (cm /sec) 

2.152x10”^ 

1.468x10”^ 

2.547x10"® 

D ^2 (om^/sec) 

3.160x10“^ 

7.795x10"^ 

7.544x10“® 

1)21 (cm /sec) 

3.483x10"^° 

8.444x1 0 "’'° 

-6 

1 .987x10 

Dgg (cm /s©c) 

2.744x10“"^ 

1 .524x10""^ 

7.766x10"® 
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The concentration dependence of the D. . ' s was also 
stronger compared to those obtained Cussler and lightfoot 
Cl965aj • In the data considering cyclohexane as solvent 
[Results ~ Set A, Table 4.4], among the main terms, ^ 
showed a large but systematic variation with concentration; 
whereas D 22 ^®niained almost constant. The cross term Dgi 
showed a trend but the variation in systematic. 

In the data considering toluene as solvent [Results - 

jfb 

Set B, Table 4.4 J, there was no discemable variation of 
the D^^’s with composition. The D 22 values were remarkably 
larger than the other terms for the first two ternary mixtures. 
Cussler and lightfoot [l965aj also reported large B 22 
values. 


The D. .'s must satisfy the criteria [Oussler and 
1 3 

lightfoot, 1965a J 


®12 ^21 ^ ° 


D^1 + B22 > ^ 


(4.2) 

(4.5) 


The first criteria was in-built in the program 
used to determine the to prevent occurrance of 

negative quantity under square root. The second criteria 
was also found to be satisfied with the data. 


; a * 


JA 
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In all the experiments, cyclohexane mass fraction 
was kept the same in the two compartments of the cell and 
thus there was no driying force for its transfer in the 
classical sense, Nevertheless, the composition changed in 
course of the runs, which indicates the non-Pickian behaviour 
as also seen from the dominance of cross terms as described 
earlier. 

One of the inherent limitation of the present technique 
is that large composition difference between the upper and 
the lower compartments is to be maintained. This is 
necessary to ensure that the change in concentration is 
much larger than the experimental errors in analysis. However, 
tha diffusivities determined thus may not represent the 
diffusivity at the arithmetic average of the compositions 
of the two Compartments due to strong non-linear nature of 
concentration dependence. One way to eliminate this error 
is to keep the composition differences very small and have 
very accurate analytical technique. The other is to 
perform experiments at different composition differences 
and extrapolate the results to zero composition difference. 
However, the effort needed is enormous in the later. 

The diffusivity data obtained in the present work 
is not adequate for obtaining a relation for the concentration 
dependence of ternary diffusivities* Development of a 
suitable correlation to account for the dependence of 
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diff usivities on concentratioii and moleculax weight of 
polymer requires more extensiye experimental work, 

4.5 Conclusions 

The ternary dif f usivit ies of polystyrene-toluene- 
oyclohexane mixture were determined using the diaphragm 
cell method. The optimum value of the diffusivity matrix 
[Dj was determined from the p ^ vs, t data employing a 
numerical search technique developed by Box. The observed 
multicomponent effects and the dependence of the diffusi- 
vities on composition were stronger than those reported ly 
Cussler and Lightfoot for the same system but differing 
only in the molecular weight of the polymer. Much more 
extensive experimental work is needed to obtain the 
dependence of diff usivities on composition and polymer 
molecular weight. 
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APPEICDIX 


A 


Table A.I : Calibration data for toluene- cyclohexane 
binary mixtures at 250 C 


""r 

W 1 

T » 

. - r 


! 

t 

! 

T 

{ 

Density, P 
(g/ml) 

0,00000 

1 .00000 


0.78269 

0.20870 

0.79130 


0.79676 

0.43303 

0.56697 


0.81442 

0.63366 

0.36634 


0.83263 

0.81163 

0.18837 


0.85059 

1 .00000 

0.00000 


0.87143 



49 


Table A. 2 : Calibration data for polystyrene-toluene- 
oyclohexane ternary mixture at 250 G 


Sample 

Ho.. 

1 

; 

f 



i 

T _ 

■""T 

* W 

{ GH 

1 . - 

1 Density, p 
(g/ml) 

- t - 

1 

0.0201 1 

0.00000 

0.97989 

0.76384 

2 

0.03971 

0.00000 

O.96O29 

0.76790 

3 

0.05984 

0.00000 

0.94016 

0.77265 

4 

0.07971 

0.00000 

0.92029 

0.77673 

5 

0.02364 

0.19593 

0.78371 

0.77684 

6 

0.04023 

0.19195 

0.76782 

0.78134 

7 

0.05989 

0.18802 

0.75208 

0.78544 

8 

0.08159 

0.18397 

0.73587 

0.78959 

9 

0.01995 

0.39202 

0.58803 

0.79137 

10 

0.03808 

0.38477 

0.57715 

0.79557 

11 

0.06016 

0.37595 

0.56394 

0.80023 

12 

0.08013 

, 0.36795 

0.55192 

0.8042 9 

13 

0.02016 

0.58790 

0.39193 

O.8O904 

14 

0.03996 

0.57602 

0.38401 

0.81349 

15 

0.06000 

0.56399 

0.37599 

0.81712 

16 

0.08019 

0.55188 

0.36792 

0.82121 

17 

0.01 989 

0.78408 

0.19602 

0.82867 

18 

0.03986 

0.76811 

0.19203 

0.83207 

19 

0.06012 

0.75191 

0.18798 

0.83614 

20 

0.07998 

0.73602 

0.18400 

0.84023 

21 

0.02025 

0.97974 

0.00000 

0,84960 

22 

0.03994 

0.96006 

0.00000 

0.86312 

23 

0.06003 

0.93997 

0.00000 

0.85699 

24 

0.07978 

0.92022 

0.00000 

0.86030 

25 

0.09007 

0.90993 

0.00000 

0.86429 

26 

0.01491 

0.98508 

0.00000 

0.85166 
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-FIG. A.2 CALIBRATION CURVE FOR POLYSTYRENE- TOLUENE-CYCLOHEXANE 
TERNARY MIXTURE AT 25°C . [Ref .Table .A.2] ^ 


Table B,1 ; Binary data for toluaie-cyclohexane mixture at 25°C 

Note ; The following data ware used for cell constant determination 
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Table B.2 : Binary data for benzene- cyclohexane mixture 

at 28°C 


Time 

(br) 

{ Compartment 
{ 1 = lower 

I U = upper 

} Density, p J 

! (g/ml) ? 

f 1 

t t 

A p P u '-P 
(g/ml) 

0 

I 

0.86512 


0 

U 

0.85016 

0.01496 

0.5 

I 

0.86245 


0.5 

u 

0.85274 

0.00971 

1.0 

I 

0.86060 


1.0 

u 

0.85487 

0.00573 

1.5 

L 

0.85973 


1.5 

u . 

0.85567 

0.00406 



Table B.3 ; Binary data for polystyrene-toluene mixture at 28°0 
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Table B.4 : Ternary raw data (Set 1) for polystyrene- 

toluene- cyclohexane mixture at 25°0 

[Ref.: Pig.B.4J 

Set 1 

Composition : Wpg = 0.»5 

Wj =* 0,90 

''oh = 


1 r 


Time 

(hr) 

f 

r 

! 

! 

f 

f 

Compartment { 

I = lower ] 

U = upper 1 

t 

Density, p J 

(g/ml) .. j 

1 

f 

^PS 

0.00 


1 

0.85477 

0.07982 

0.00 


U 

0.84369 

0.02002 

2.00 


L 

0.85351 

0.07556 

2.00 


U 

0.84490 

0.02606 

5.25 


1 

0.85272 

0.07180 

3.25 


U 

0.84502 

0.02882 

4.50 


L 

0.85297 

0.07223 

4.50 


U 

0.84540 

0.02963 

5.75 


I 

0.85268 

0.06892 

5.75 


u 

0.84581 

0.03326 
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Table B.5 : Temar 7 raw data (Set 2) for polystyrene- 

toluene-cyclohexane mixture at 25 °C 

CEef. : Pig. B.5 J 


Set 2 


Composition ; 


Wpg = 0.05 

Wj = 0.70 

WqH = 0*^5 


T r 


Time 

(hr) 

{ Compartment 

1 1 = lower 

{ U = upper 

- f - - - 

} Density, P | 

; (g/ml ) ; 

1 ; 

1 1 

^PS 

0 

I 

0.83314 

0,07932 

0 

U 

0.82297 

0.01917 

1 

1 

0.83299 

0.07699 

1 

U 

0.82509 

0.02194 

2 

L 

0.83276 

0.07560 

2 

U 

0.82289 

0.02360 

3 

L 

0.83268 

0.07323 

3 

U 

0.82356 

0.02536 

4 

1 

0.83221 

0.07319 

4 

U 

0.82394 

0.02685 

5 

1 

0.83218 

0.07094 

5 

u 

0.82391 

0,02847 
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Table B,6 ; Ternary raw data (Set 3) for polystyrene- 

toluene-cyclohexane mixture at 25 °0 

[Ref , : Pig,B.6 ] 

Set 3 


Composition : 


Wpg = 0.05 
Wp = 0.475 
^GH ~ ^*475 


Time 

(hr) 



Compartment 

1 L = lower 

1 U = upper 

..-p,,..— 

* Density, p ’ 

(g/ml ) 1 

t f 

..I »' 


0 

1 

0.81119 

0.07897 

0 

U 

0.801 42 

0.01960 

1 

1 

0.81093 

0.07837 

1 

U 

0.80180 

0.02192 

2 

I 

0.81019 

0.07635 

2 

u 

0.80221 

0.02397 

3 

1 

0.81031 

0.07465 

3 

u 

0.80224 

0.02518 

4 

I 

0.80960 

0.07282 

4 

u 

0.80329 

0.02678 

5 

1 

0,80937 

0.07250 

5 

u 

0.80334 

0.0281 1 
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Tatle B.g : Smoothanad ternary data (Set 3) for polystyraae-toluane-cyclohexane mixture 
at 25*0 

Note; The following data are read from graphs in lig,B»6 

Set 3 

Composition : Wpc, = 0.05 
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Table B«11 t Ternary P . vs, t data (Set 2) for polystyrene-toluaie-cyclohexane mixture 
at 25°C. 

Note : The following data were calculated using Table B.8 

Set 2 

Composition: Wp^ = 0,05 
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10800 U 0,0208905 0. 5941484 0.2085812 0.0400686 -0.0267928 -0.0044158 

14400 L 0.0597898 0.5691984 0.2033419 

14400 U 0.0221181 0.5933454 0.2084265 0.0376717 -0.0241470 -0.0050847 

18000 L 0.0590417 0.5702535 0.2029249 

18000 U 0.0234611 0.5917241 0.2089148 0.0355805 -0.0214706 -0.0059899 


B.12 : Ternary P j_ vs, t data (Set 3) for polystyrene-toluene-oyclohexane mixture 
at 25°C 

Note ; The following data were calculated using Table B,9 
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FIG B1 BINARY DATA FOR TOLUENE - CYCLOHEXANE 
MIXTURE AT 28 »C FOR DIAPHRAGM 1 
[Ref. -.Table B.l] 





FIG B-2 BINARY DATA FOR TOLUENE - CYCLOHEX ANE 
MIXTURE AT 28°C FOR DIAPHRAGM 2 . 

[Ref. -.Table B.l] 



Run-time {hr ) 

FIG B.3 BINARY DATA FOR BENZENE - CVLOHEXANE 
MIXTURE AT ZS^C [Ref.: Table B.z] 




Mass fr. polys tyrene^Wps Mixture density 
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-CYCLOHEXANE MIXTURE (SET 1 ) AT 25»C 
[Ref.; Table B.4] 




Mass fr. polys trene^^Wps Mixture density 
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R6. 0-5 



TERNARY DATA FOR POLYSTYRENE- TOLUENE 
[Ref.;Tol>le 0.5] ' 





Mass fr. polystyrene ,Wps Mixture density (g/ml ) 




FIG. e.6 TERNARY DATA FOR POLYSTYRENE -TOLUENE 
-CYCLOHEXANE MIXTURE (SET 3) AT aS^C 
[Rcf.--Table B.6] 




